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ABSTRACT 

This  thesis  presents  a  comparison  of  three  Magnetic 
Anomaly  Detection  (MAD)  models:  a  cross-correlation 
detection  model,  a  square  law  detection  model,  and  a  model 
referred  to  as  the  OPTEVFOR  detection  model.  FORTRAN  and 
BASIC  programs  for  the  three  detection  models  are  included 
in  this  thesis.  The  programs  yield  detection  probabilities 
for  straight  line  encounters.  Magnetic  signal  values  for 
the  straight  line  encounters  are  an  additional  output. 
Plots  of  lateral  range  curves  and  magnetic  signal  values  are 
presented.  A  discussion  of  the  required  parameters  is 
included  in  the  thesis  to  facilitate  the  use  of  the 
programs.  The  parameters  that  were  considered  in  the 
comparison  of  the  three  detection  models  are:  magnetic 
noise,  aircraft  and  submarine  headings,  submarine 
displacement,  and  the  vertical  separation  between  submarine 
and  aircraft. 
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I.  INTRODUCTION 


This  thesis  presents  a  comparison  of  three  Magnetic 
Anomaly  Detection  (MAD)  models.  The  comparison  is  in  terms 
of  probabilities  of  detection  that  were  computed  using  the 
models.  Two  of  the  models,  the  cross-correlation  model  and 
the  square  law  model,  have  been  used  to  model  sonar 
detection  [Ref.  1:  pp.  343-357  J.  The  third  model,  referred 
to  as  the  OPTEVFOR  model,  is  a  slant  range  threshold 
detection  model.  The  results  of  the  comparisons  are 
presented  in  graphical  and  tabular  form.  In  addition,  plotB 
of  magnetic  signals  for  selected  lateral  ranges  and  noise 
levels  are  shown.  The  effects  of  noise,  aircraft  and 
submarine  headings,  submarine  displacement,  and  vertical 
separation  are  also  indicated. 

The  models  were  implemented  using  the  fortran  and  BASIC 
programs1  that  are  listed  in  Appendix  A.  For  those 
interested  in  using  the  programs  for  other  investigations, 
an  input  parameter  discussion  is  provided  in  Chapter  3.  To 
use  the  FORTRAN  program,  the  user  specifies  the  input 
parameters  in  an  input  file.  After  execution  of  the 
program,  an  output  file  is  generated  that  contains 

^he  programs  are  based  on  an  unpublished  BASIC  program 
by  R.N.  Forrest  for  an  H.P.-  05  microcomputer. 


probabilities  of  detection  for  each  of  the  three  models.  In 
addition,  magnetic  signal  values  and  magnetic  signal  values 
plus  random  magnetic  noise  values  for  one  of  the  encounters 
generated  by  the  program  are  included  in  the  output  file. 
An  IBM  GRAFSTAT  graphical  package  was  used  to  produce  the 
graphics  in  this  thesis. 

To  use  the  BASIC  program,  the  user  must  interactively 
enter  the  input  parameter  values  for  each  encounter.  After 
execution  of  the  program,  an  optional  hardcopy  printout 
supplies  the  input  parameter  values  and  a  table  of  detection 
probabilities  for  each  of  the  three  models  (see  Appendix  A). 
Following  this,  lateral  range  curves  are  displayed  to  the 
user  for  immediate  observation.  A  typical  program  run 
producing  21  detection  probabilities  for  each  model  requires 
approximately  10  minutes  of  computing  time  on  an  Atari  800 
microcomputer. 


the  limiting  detection  capability  for  an  automatic  system 
that  does  not  use  information  about  the  shape  of  the 
magnetic  signal. 


B.  OPTEVFOR  MAD  DETECTION  MODEL 

The  OPTEVFOR  model  is  described  by  Forrest  [Ref.  3:  pp. 
7-8].  In  characterizing  the  submarine  magnetic  signal  as  a 
simple  dipole  signal,  the  U.S.  National  Defense  Research 
Committee,  [Ref.  4:  p.  20],  reports  that  the  magnetic  signal 
of  the  submarine  "varies  as  the  inverse  cube  of  the  distance 
from  the  source".  In  an  OPTEVFOR  report  [Ref.  5:  p.  1  ,  end. 
1],  the  results  of  a  regression  analysis  on  empirical  peak 
to  peak  signal  output  against  slant  range  between  submarines 
and  aircraft  are  reported.  These  results  also  suggested 
this  inverse  cube  relationship  for  the  magnetic  signal. 
This  relationship  is  the  basis  for  the  OPTEVFOR  detection 
model . 

The  model  has  a  deterministic  mode  and  a  stochastic 
mode,  each  of  which  involves  the  following  parameters:  the 
submarine  magnetic  moment  (M),  an  Operator  Recognition 
Factor  (ORF),  the  average  peak  to  peak  magnetic  noise  (N)  in 
the  operating  area,  and  a  slant  range  (R).  The  relationship 
between  these  quantities  is  given  by: 


C  M 

(ORF)  N 


1/3 


( eqn  2.1) 


The  value  of  the  constant  c  is  0.^0  for  M  in  oersted 
centimeters3,  R  in  meters,  and  N  in  gamma. 

In  the  deterministic  mode,  detection  occurs  and  only 
if  the  aircraft's  slant  range  from  the  submarin  at  CPA  is 
less  than  or  equal  to  R.  This  mode  yields  a  rectangular 
("cookie  cutter")  lateral  range  curve  with  the  probability 
of  detection  equal  to  1  for  an  encounter  where  the  slant 
range  at  CPA  is  less  than  or  equal  to  R,  and  0  whon  it  is 
greater  than  R. 

The  stochastic  mode  allows  a  more  uncertain  approach  to 
detection  by  allowing  a  gradual  rise  in  probability  of 
detection  as  the  slant  range  at  CPA  decreases.  In  this  mode 
one  sets  the  probability  of  detection  at  R  equal  to  50 
percent,  and  the  lateral  range  curve  is  given  by  pd  »  <|)(x); 
where  it  is  understood  that  ^  is  the  standard  normal 
cumulative  distribution  function  and  x  is  determined  by  the 
following  equation: 

R  -  CPA 

X  a  - 

(AL)  R  .  (eqn  2.2) 
In  this  equation,  CPA  is  the  magnitude  of  the  slant  range 
distance  at  CPA,  and  R  is  the  calculated  range  from  Equation 
2.1.  The  product  { AL)R  represents  a  standard  deviation. 
The  value  of  AL  can  be  considered  to  be  determined  by  "the 
combined  uncertainty  and  variability  in  the  values  of  Mr  N, 
and  ORF"  [Ref.  3:  p.  8].  Two  values  of  AL  (.20  and  .01)  are 


shown  in  Figure  2.1.  If  empirical  data  was  available,  the 
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value  of  AL  could  be  chosen  to  provide  a  best  fit  to  the 
observed  res-’ Its.  Note,  as  AL  approaches  0,  the  stochastic 
mode  approximates  the  deterministic  mode. 


AL  ■  0.80  AL  -  0,01 


Figure  2.1.  Lateral  Range  Curves  for  Different  Values 
of  AL. 


III.  INPUT  PARAMETERS 


The  input  parameters  for  the  FORTRAN  program  are  all 
contained  in  one  input  file.  This  allows  parameter  values 
to  be  easily  changed  without  recompiling  the  main  program  or 
subroutines.  Also,  with  a  few  changes,  this  program  could 
be  altered  to  operate  in  conju*  ion  with  a  larger  program 
to  yield  a  probability  of  detection  on  an  individual  MAD 
run. 

The  input  parameters  are  divided  into  four  areas  for 
discussion.  They  ares  (1)  sample  interval,  (2)  earth 
magnetic  field,  (3)  submarine  moments,  and  (4)  other 
inputs. 

A.  SAMPLE  INTERVAL 

The  choice  of  a  sample  interval  is  discussed  by  Forrest 
[Ref.  2:  pp.  27-30  ].  In  the  program,  the  total  observation 
time  in  seconds  over  which  the  samples  are  taken  is  entered 
in  T7.  This  time  should  be  long  enough  to  encompass  a 
"complete  signal"  at  the  maximum  expected  detection  slant 
range . 

As  the  slant  range  from  the  submarine  to  the 
magnetometer  increases,  the  distance  over  which  a 
significant  magnetic  signal  is  present  at  the  magnetometer 
also  increases.  Figure  3.1  graphically  shows  the  difference 
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in  the  amount  of  time  that  a  signal  is  present  for  slant 
ranges  of  200  meters  and  805  meters.  In  this  thesis,  the 
total  time  for  a  straight  line  encounter  is  assumed  to  be  20 
seconds.  As  can  be  be  seen  from  Figure  3.1  ,  a  20  second 
interval  adequately  covers  the  significant  portion  of  the 
magnetic  signal  for  an  805  mmeter  slant  range  at  CPA. 


MO  MC1IRS  CM  SUNT  MMOC  805  METERS  CM  SLANT  RANGE 


SECONDS  SECONDS 

Figure  3.1.  Magnetic  Signals  for  Slant  Ranges  of  200  Meters 
and  805  Meters. 

The  time  between  samples  is  set  equal  to  the  reciprocal 
of  twice  the  upper  bandpass  filter  frequency  of  the  MAD 
sensor.  A  value  of  0.9  Hz  was  suggested  for  use  by  Texas 
Instruments  [Ref.  6:  p.  112]  as  an  upper  bandpass  filter 
limit  in  a  discussion  on  the  effects  of  noise  on  a  MAD 
system.  This  value  yields  a  time  interval  between  samples 
of  0.55  seconds. 

The  sample  interval  length  and  the  false  alarm  rate  (the 
expected  number  of  false  alarms  per  hour)  determine  the 
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false  alarm  probability.  The  false  alarm  rate  (F2)  is 
assigned  a  value  of  3  based  on  a  report  by  OPTEVFOR  [Ref.  5: 
p.  2.1]. 

B.  EARTH  MAGNETIC  FIELD 

Input  values  for  the  earth  magnetic  field  intensity  and 
inclination,  or  dip  angle,  may  be  taken  from  two  Defense 
Mapping  Agency  Hydrographic  Center  charts,  [Refs.  7  and  8 
respectively],  or  approximated  by  using  a  program.  If  chart 
values  are  entered,  the  earth  field  intensity  must  be  in 
units  of  gamma  and  the  inclination  in  decimal  degrees.  The 
program  used  to  determine  the  intensity  of  the  earth  field 
and  inclination  is  based  on  a  simple  dipole  field  model  that 
ii  described  by  Forrest  [Ref.  9:  pp.  39-43]. 

Table  III-1  displays  the  program  output  values  of 
inclination  in  decimal  degrees  and  earth  magnetic  field  in 
gamma  for  selected  geographic  locations.  In  addition, 
corresponding  values  obtained  from  the  Defense  Mapping 
Agency  Hydrographic  Center  Charts  Number  30  and  Number  39 
are  also  displayed.  The  last  three  columns  are  the  average 
slant  range  in  meters  at  which  a  50  percent  probability  of 


detection  is  obtained  for  the  three  program  detection 
models.  The  program  input  parameters  for  these  slant  ranges 
were  the  same  as  the  base  case,  except  for  the  following 
differences:  a  sample  interval  time  of  40  seconds,  aircraft 
and  submarine  headings  of  0  degrees,  and  a  submarine 


displacement  of  7,000  tons.  The  correlation  between  the 
slant  ranges,  comparing  the  chart  values  and  model  values, 
was  found  to  be  95  to  96  percent  for  the  three  models.  This 
suggests  that,  even  though  differences  exist  between  the 
chart  values  and  model  values,  there  is  a  high  degree  of 
correlation  in  the  final  output. 

A  limitation  to  the  simple  dipole  field  model  is  that 
it  does  not  give  an  angle  of  declination  (variation)  with 
sufficient  accuracy.2  As  a  result,  all  headings  entered 
into  this  program  must  be  in  magnetic  degrees.  The  Phoenix 
Corporation  [Ref.  10:  pp.  24-25)  reports  on  geomagnetic 
field  models  that  can  represent  the  earth  field  "with 
overall  accuracies  better  than  approximately  150-200  gammas 
in  magnitude  and  .2*  in  direction  of  the  field."  This 
degree  of  accuracy  is  not  needed  for  this  program,  but  a 
simplified  version  of  one  of  these  models  that  provided 
satisfactory  angles  of  declination  would  be  beneficial  if 
the  program  were  to  be  incorporated  into  a  larger  model  that 
utilized  true  headings  as  inputs. 

C.  SUBMARINE  MAGNETIC  DIPOLE  MOMENT 

If  a  submarine’s  magnetic  dipole  moment  is  known  for  the 
geographical  location  and  the  submarine's  magnetic  heading, 

2Private  communication  from  R.N.  Forrest,  who 
investigated  the  use  of  the  simple  dipole  model  for  this 
purpose. 


the  following  values  may  be  entered  in  the  program:  (1)  P, 
its  magnitude  in  oersted  centimeters  cubed,  (2)  A,  its 
direction  in  decimal  degrees  relative  to  magnetic  north,  and 
(3)  B,  its  depression  angle  from  the  horizontal  in  decimal 
degrees.  If  it  is  not  known,  these  values  must  be 
calculated  for  a  specific  location  and  magnetic  heading.  A 
program  is  included  in  the  main  program  that  can  be  used  to 
calculate  these  values.  The  program  is  based  on  a  model 
described  by  Forrest  [Ref.  9:  pp.  35-38  ].  The  input  to  the 
program  is  submarine  displacement  in  tons.  The  program  also 
contains  coefficients  which  relate  displacement  to  magnetic 
moment.  The  values  used  in  the  program  are  based  on  values 
cited  by  Texas  Instruments  [Ref.  6:  p.  4]. 

The  past  history  of  the  submarine  is  represented  by  the 
permanent  longitudinal,  transverse,  and  vertical  moments  of 
the  submarine  (M4,  M5  and  M6  in  oersted  centimeters  cubed). 
For  the  examples  in  this  thesis,  it  was  assumed  that 
effective  deperming  had  been  performed  and  program  default 
values  of  zero  were  used. 

D.  OTHER  PARAMETERS 

1  .  Headings  and  Speeds 


Since  the  simple  dipole  earth  field  model  used  by 
the  program  does  not  produce  accurate  angles  of  declination, 
magnetic  headings  are  required.  In  addition,  the  headings 


must  be  in  decimal  degrees.  The  input  parameters  for 
submarine  speed  and  aircraft  speed  are  entered  in  knots. 

2.  Noise 

The  magnetic  noise  is  assumed  to  be  such  that 

adjacent  magnetic  noise  samples  are  independent.  This 

assumption  is  based  in  part  on  the  filtering  that  is 

performed  on  the  magnetic  signal  by  the  processing  system  in 

a  MAD  detection  sensor.  The  standard  deviation  of  the  noise 

in  units  of  gamma  is  the  value  entered  into  SI.  This  value 

can  be  approximated  from  operational  data  by  talcing  from 

one-fourth  to  one-sixth  of  the  measured  peak  to  peak 
) 

magnetic  noise.  [Ref.  2:  pp.  28-29} 

The  OPTEVFOR  detection  model  incorporates  a  value  of 
average  peak  to  peak  magnetic  noise  (N)  in  the  inverse  cube 
law  calculation.  In  the  program,  the  value  of  N  is 
determined  by  multiplying  the  SI  entry  by  four. 

3.  Operator  Recognition  Factor  (ORF) 

The  ORF  is  the  value  of  the  ratio  of  magnetic  signal 
to  magnetic  noise  for  which  the  average  operator  would 
detect  a  signal  50%  of  the  time  in  the  presence  of 
background  noise  for  a  false  alarm  rate  of  3  per  hour.  An 
ORF  value  of  3  was  suggested  for  use  by  OPTEVFOR  [Ref.  5:  p. 
4.12]. 

4.  Distance  Parameters 


Two  parameters,  R8  and  N7,  are  used  to  define  the 
points  plotted  on  the  lateral  range  curves.  R8  is  the 


maximum  positive  value  of  the  lateral  range  in  meters  for 
which  a  lateral  range  curve  value  is  to  be  computed.  N7 
represents  the  number  of  lateral  range  curve  values  that  are 
to  be  computed  from  the  maximum  lateral  range  to  zero 

lateral  range. 

The  vertical  separation  (Z)  is  the  sum  of  the  submarine 
depth  and  aircraft  altitude  in  meters. 


IV.  RESULTS 


Program  outputs  of  the  three  models  for  a  set  of  base 
case  conditions  are  presented  in  this  section.  Outputs  for 
variations  from  the  base  case  are  also  presented.  The 
lateral  range  of  an  encounter  (the  horizontal  separation 
between  the  submarine  and  magnetometer  when  the  magnetometer 
is  at  CPA)  for  a  50%  probability  of  detection  is  used  as  a 
measure  of  comparison.  Signal  and  signal  plus  "noise" 
traces  for  several  cases  are  presented.  The  traces  are 
based  on  the  signal  and  noise  models  that  are  part  of  the 
cross-correlation  and  square  law  models.  These  idealized 
signal  traces  appear  to  have  the  characteristics  of  actual 
signal  traces.  This  suggests  that  the  signal  and  noise 
models  might  be  used  for  training  purposes. 

A.  BASE  CASE 

The  base  case  conditions  are  listed  in  Table  IV- 1.  The 
table  is  ordered  in  the  same  manner  that  the  values  are  read 
into  the  program.  An  annotation  of  each  entry  is  included 
for  clarity. 

Figure  4-1  presents  the  lateral  range  curves  for  the 
base  case.  Points  on  the  lateral  range  curves  are  indicated 
by  the  first  letter  of  the  name  of  the  model  from  which  they 
ware  derived.  The  slight  asymmetry  of  the  cross-correlation 
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detection  model  and  square  law  detection  model  curves  is 
reflective  of  the  shape  of  the  signals  that  are  'received' 
in  these  models. 

Table  IV- 1 .  Input  Parameters  for  the  Base  Case 


1  .8 
20.0 
3.0 
0 

30.0 

60.0 

45.0 

10.0 

315.0 

220.0 

0 

0 

0 

4000.0 

200.0 

0.1 

1500.0 

50.0 

3.0 

0.2 

0025 


twice  the  upper  bandpass  limit  in  seconds 
sampling  time  interval  in  seconds 
false  alarms  per  hour 
Enter  inclination  (1  ■  yes,  0  ■  no)? 
area  of  operation  latitude  in  decimal 
degrees 

area  of  operation  longitude  in  decimal 

degrees 

submarine  magnetic  heading  in  decimal 
degrees 

submarine  speed  in  knots 

aircraft  magnetic  heading  in  decimal 

degrees 

aircraft  speed  in  knots 

Enter  submarine  moment  (1  ■  yes,  0  -  no)? 

Enter  earth  field  (1  ■  yes,  0  ■  no)? 

Enter  submarine  perm  moments  (1«  yes, 

0  ■  no)? 

submarine  displacement 
vertical  separation  in  meters 
noise  (standard  deviation)  in  gamma 
maximum  lateral  rango  in  meters 
divisions  of  lateral  range 
ORE  (Operator  Recognition  Factor) 
variability  factor  for  OPTEVFOR  model 
lateral  range  iteration  number  for  the 
magnetic  signal  and  signal  plus  noise  in 
the  output  file 


Table  IV-2  lists  lateral  detection  ranges  and 
corresponding  slant  detection  ranges  at  CPA  for  a 
probability  of  detection  equal  to  50  percent  for  the  cross¬ 
correlation  and  square  law  detection  models.  An  equivalent 
ORF.  value  for  each  model  is  also  listed.  Due  to  the 
asymmetry  of  the  lateral  range  curves  for  the  cross- 
correlation  and  square  law  models,  the  average  of  the  two  50 


-1000 


0 

METERS 


Figure  4-1.  Lateral  Range  Curves  of  the  Cross-Correlation 
(C),  Square  Law  (S),  and  OPTEVFOR  (0)  Models  for  the  Base 
Case. 


Table  IV-2.  The  Lateral  Detection  Ranges,  Slant  Detection 
Ranges,  and  ORF's  of  the  Three  Models  for  the  Base  Case. 

Lateral  Detection  Slant  Detection  ORF 
Range  (meters)  Range  (meters) 

Cross-Correlation  885  907  .21 

Square  Law  685  714  .44 

OPTEVFOR  318  376  3 

percent  detection  ranges  was  used  as  the  lateral  detection 
range.  The  equivalent  ORF  values  for  the  cross-correlation 


and  square  law  detection  models  were  calculated  using  the 
slant  detection  range  values  with  the  following  equation, 
which  was  obtained  from  Equation  2.1: 

c  M 

ORF  »  ---- 

R3  N  .  eqn  4.1 

For  the  base  case,  the  magnitude  of  the  submarine  field  (M) 
at  the  submarine  is  6.35  x  108  orested  cm3,  the  noise  (N)  is 
.4  gamma,  and  the  value  of  the  constant  (c)  is  .1.  This 
suggests  that,  in  order  to  detect  a  magnetic  signal  50 
percent  of  the  time  with  a  false  alarm  rate  of  3  per  hour, 
the  magnetic  signal  to  magnetic  noise  ratio  should  be  .21 
for  an  ideal  cross-correlation  detector  and  .44  for  an  ideal 
square  law  detector. 

Using  the  ORF  values,  the  cross-correlation  and  square 
law  detection  models  can  be  used  to  describe  the  performance 
of  an  operator.  To  do  this,  a  modified  value  of  the 
standard  deviation  (o')  of  the  input  noise  can  be  used.  The 
modified  value  is  oqual  to  (ORF)  (o')/. 21  for  the  cross¬ 
correlation  detection  model  and  (ORF) ( «■*)/. 4 4  for  the  square 
law  detection  model.  With  these  modifications,  the  two 
models  can  be  used  to  describe  the  detection  capability  of 
an  operator  with  a  specified  ORF.  An  example  of  a  lateral 
range  curve  with  the  modified  noise  standard  deviation  for 
an  ORF  of  3  is  presented  in  Figure  4.2  for  each  model. 
These  curves  are  comparable  to  the  lateral  range  curve  for 
the  OPTEVFOR  model  that  is  shown  in  Figure  4.1. 
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Figure  4.2.  Cross-Correlation  and  Square  Law  Lateral  Range 
Curves  to  Describe  the  Performance  of  an  Operator  with  an 
ORF  of  3. 
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The  automatic  MAD  system  manufactured  by  Canada's  CAE 
Electronics  Ltd.  is  expected  to  produce  a  50  percent  in¬ 
crease  in  detection  slant  range  [Ref.  11].  Using  the 
detection  slant  range  for  the  OPTEVFOR  model  of  376  meters, 
a  50  percent  improvement  would  yield  a  detection  slant  range 
of  564  meters.  The  ORF  for  a  detection  system  with  this 
capability  would  be  .88.  The  cross-correlation  and  the 
square  law  detection  models  could  be  used  to  yield  lateral 
range  curves  for  a  system  with  an  ORF  of  .88  by  using  a 
noise  standard  deviation  equal  to  .88  (O')/. 21  and  .88  (o**)/ 
.44  respectively.  Figure  4.3  shows  the  lateral  range  curves 
of  the  two  detection  models  with  a  50  percent  improvement  in 
slant  range  detection.  Note,  with  the  modified  noiae 
standard  deviations,  the  models  are  essentially  equivalent 
for  the  cases  considered. 
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Figure  4.3.  The  Cross -Correlation  and  Square  Law  Models  to 
Describe  LRC's  for  the  CAE  Automatic  Detection  System. 


SK3NAL  AT  CPA  or  0  MITERS 


SIGNAL  +  MOST  AT  CPA  OT  0  MOER9 


Figure  4„4.  Magnetic  Signal  and  Magnetic  Signal  Plus 
Magnetic  Noise  at  a  Lateral  Range  at  CPA  of  0  Maters  for  the 
Base  Case. 


Figures  4.4  and  4.5  present  the  magnetic  signal  and  a 
representation  of  magnetic  signal  plus  magnetic  noise  that 
would  be  received  under  the  base  conditions  by  a 
magnetometer  with  a  lateral  range  of  0  meters  and  of  780 
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meters.  The  signal  plus  noise  trace  was  generated  from 
signal  plus  noise  values  obtained  by  adding  a  signal  value 
to  a  gaussian  noise  value.  The  gaussian  noise  value  was 
generated  by  multiplying  the  standard  deviation  of  the  input 
noise  by  a  pseudo  normal  random  number  from  a  population 
with  mean  0  and  variance  1.  The  pseudo  normal  random  numbers 
were  generated  using  LLRANDOMII,  a  resident  program  at  the 
Naval  Postgraduate  School  computer  [Ref.  12:  p.  2.2]. 


StONM.  AT  CCA  Of  TOO  METERS  SGNAL  +  NOSE  AT  CCA  Of  780  METERS 


SECONDS  SECONDS 


Figure  4.5.  Magnetic  Signal  and  Magnetic  Signal  Plus 
Magnetic  Noise  at  a  Lateral  Range  at  CPA  of  780  Meters  for 
the  Base  Case. 

The  magnitude  of  the  magnetic  signal  shown  in  Figure  4.4 
is  very  large  in  comparison  to  the  background  noise.  The 
peak  to  peak  signal  to  noise  ratio  is  approximately  14  to  1. 
An  operator  would  have  little  difficulty  identifying  the 
signa.u  in  this  signal  plus  noise  trace. 

Conversely,  the  magnetic  signal  shown  in  Figure  4.5  is 
small  compared  to  the  background  noise.  The  peak  to  peak 


detection  for  the  lateral  range,  of  780  meters  are:  .95  for 
the  cross-correlation  detection  model,  .28  for  the  square 
law  detection  model,  and  0  for  the  OPTEVFOR  detection  model. 
It  seems  apparent  that  an  operator  would  have  a  difficult, 
if  not  impossible,  time  in  detecting  this  signal  at  a 
reasonable  false  alarm  rate. 

B.  DIFFERENT  NOISE  INPUTS 

The  first  variation  on  the  base  case  shows  the  effect  of 
different  noise  inputs.  The  standard  deviation  (O')  of  the 
peak  to  peak  noise  is  the  input  parameter  that  is  varied. 
Table  IV-3  lists  the  different  O'  values  and  the 
corresponding  lateral  detection  ranges. 

Table  IV-3.  The  Effect  of  Noise  on  Detection  Range. 

Standard  Deviation  Lateral  Detection  Range  in  Meters 
of  Noise  in 

Gamma  Cross - 


Correlation 

Square  Law 

OPTEVFOR 

.005 

2  2  50 (  2  2  5  9)* 

1  7  9  2 (  1  803)* 

1000  (1020)* 

.01 

1832  (1843) 

1440 

(1460) 

782  (807) 

.05 

1110  (1128) 

863 

(890) 

427  (472) 

.1 

885  (907) 

685 

(714) 

318  (375) 

.5 

512  (550) 

382 

(431  ) 

90  (219) 

*The 

numbers  in  parentheses 

are  the  slant  range 

distances  in  meters.  The  vertical  separation  is  200  meters. 

Figure  4^6  displays  lateral  range  curves  for  the  three 
models  when  the  standard  deviation  of  the  noise  is  .01 
gamma.  These  three  curves  show  an  increase  in  lateral 


detection  range  over  the  base  case.  Note  that  the  asymmetry 
of  the  cross-correlation  and  square  law  detection  models  is 
more  apparent  in  Figure  4.6  than  it  was  in  Figure  4.1. 

Figure  4.7  displays  the  magnetic  signal  (which  is  the 
same  as  the  signal  in  Figure  4.5)  and  the  magnetic  signal 
plus  magnetic  noise  at  a  horizontal  distance  of  780  meters 
when  the  magnetometer  is  at  CPA.  The  signal  to  noise  ratio 
is  3.5.  The  figure  suggests  that  a  MAD  operator,  in  this 
case,  should  have  the  ability  to  detect  a  signal  at  780 
meters  lateral  range  with  a  satisfactory  false  alarm  rate. 


Figure  4.C.  Lateral  Range  Curves  for  the  Three  Models  with 
the  Standard  Deviation  of  the  Noise  Set  to  .01  Gamma. 
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SKm.  4  NOISE  AT  CPA  OF  780  METERS 


Figure  4.7.  Magnetic  Signal  and  Magnetic  Signal  Plus 
Magnetic  Noise  with  the  Standard  Deviation  of  Noise  »  .01 
Gamma  at  780  Meters  Lateral  Range. 

C.  DIFFERENT  HEADINGS 

The  headings  of  a  submarine  and  an  aircraft  in  an 
encounter  have  an  effect  on  detection  ranges.  The  effect  of 
different  headings  was  investigated  using  the  square  law 
detection  model,  and  the  results  in  terms  of  lateral 
detection  ranges  are  presented  in  Table  IV-4.  This  table 
suggests  that  a  submarine  should  choose  a  magnetic  heading 
of  either  East  or  West,  and,  for  an  encounter,  an  aircraft 
should  also  choose  a  magnetic  heading  of  East  or  West. 

Figure  4.8  shows  lateral  range  curves  for  a  submarine 
heading  North  and  an  aircraft  heading  East.  In  this  case, 
both  the  cross-correlation  and  square  law  detection  model 
lateral  range  curves  display  noticeable  asymmetry.  The 
OPTEVFOR  detection  model  lateral  range  curve  is  symmetric 
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Table  IV- 4 .  Square  Law  Lateral  Detection  Ranges  for 
Different  Submarine  and  Aircraft  Magnetic  Headings 

Aircraft  Submarine  Headings  (magnetic) 

Headings 

(magnetic) 

0  45  90  135  180  225  270  315 


0  700  650  498  640  686  640  498  650 

45  724  637  505  670  713  624  501  676 

90  730  654  527  654  730  646  524  646 

135  712  685  505  636  724  682  502  624 

180  685  640  498  650  700  650  498  640 

225  712  624  501  672  724  636  506  672 

270  730  646  51  9  646  730  654  533  654 

315  724  685  501  624  713  682  505  636 


Figure  4.8.  Lateral  Range  Curves  for  the  Submarine  Heading 
North  and  the  Aircraft  Heading  East. 
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but,  like  the  curves  for  the  other  models,  it  shows  an 
increase  in  detection  ranges  over  those  for  the  base  case 
(where  the  submarine  is  heading  NE  and  the  aircraft  is 
heading  NW ) . 

The  APAIR  MOD  2.6  (Ref.  13:  p.  83]  simulation  uses  a  MAD 
detection  model  that  accounts  for  the  change  in  a 
submarine's  magnetic  moment  (which  is  dependent  on  changes 
in  submarine  heading)  by  using  a  parameter  labeled  DFACTR 
(degradation  factor  for  heading).  In  the  model,  D  (a 
modified  slant  range  at  CPA)  determines  the  probability  of 
detection.  The  value  of  D  is  determined  using  the  following 
relation: 

D  -  DC  (1  -  DFACTR  x  A),  eqn.  4.1 

where  DC  is  the  slant  range  at  CPA  and  A  is  the  acute  angle 
in  decimal  uegrecs  between  the  submarine  heading  and  an 
East-West  bearing.  The  probability  of  MAD  detection  is 
determined  from  a  table  of  probability  of  detection  against 
slant  range.  A  uniform  (0,  1)  random  number  is  drawn  to 
determine  whether  or  not  the  submarine  is  detected.  The 
average  slant  detection  ranges  (computed  from  Table  IV-4 , 
where  the  vertical  separation  is  200  meters)  for  submarine 
headings  of  North  and  East  are  741  meters  and  545  meters 
respectively.  These  ranges  yield  a  value  of  ,003  for 
DFACTR.  The  average  slant  detection  range  from  Table  IV-4 
for  a  submarine  heading  of  NE  is  68  2  meters;  however,  the 
slant  range  determined  by  a  modified  slant  range  of  545 


meters  and  a  DFACTR  «  .003  is  643  motors.  If  sin  A  instead 
of  A  is  used  in  Equation  4.1,  then  DFACTR  is  .265  and  the 
slant  detection  for  a  submarine  heading  NE  is  670  meters. 
Since  this  is  only  a  single  data  point  and  there  is  no 
supporting  operational  data,  the  modification  is  not 
proposed  as  one  that  should  be  adopted.  However,  this 
cursory  analysis  does  indicate  a  way  in  which  the  programs 
presented  in  this  thesis  might  be  used  by  others. 

Table  IV-5  lists  lateral  ranges  for  P(det)  equal  to  50 
percent  for  3  submarine/aircraft  heading  combinations.  The 
cross-correlation  and  OPTEVFOR  detection  model  results  show 
the  same  relationship  as  the  results  of  the  square  law 
detection  model. 

Table  IV-5.  Lateral  Ranges  for  P(det)  -  .50  in  MeterB  for 
the  Three  Detection  Models. 


Submarine 

45 

0 

90 

Aircraft 

315 

90 

0 

Cross-Correlation 

885 

934 

754 

Square  Law 

685 

730 

498 

OPTEVFOR 

318 

358 

230 

For  the  detection  ranges  reported  by  OPTEVFOR  [Ref.  5: 
p.  5.1],  the  effect  of  different  headings  was  averaged  out. 
That  is,  measurements  were  taken  from  the  16  possible 
combinations  of  the  4  cardinal  submarine  and  aircraft 
headings  in  equal  numbers  and  then  averaged  to  yield  an 
average  slant  detection  range.  But,  as  shown  in  Tables  IV-4 
and  IV-5,  the  models  show  significant  variability  in  lateral 
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detection  range  for  different  submarine  and/or  aircraft 
headings. 

Figure  4.9  is  included  to  show  the  lateral  range  curves 
when  the  submarine  is  headed  East  and  the  aircraft  is  headed 
North.  These  lateral  range  curves  give  the  minimum  lateral 
detection  ranges  for  the  different  heading  combinations. 
Also,  for  the  cross-correlation  and  square  law  detection 
models,  the  lateral  range  curves  are  fairly  symmetric. 


Figure  4.9.  Lateral  Range  Curves  for  the  Submarine  Headed 
East  and  the  Aircraft  Headed  North. 
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D.  SUBMARINE  DISPLACEMENT 


The  submarine  magnetic  dipole  moment  program  within  the 
main  program  is  used  to  calculate  a  submarine's  induced 
magnetic  moments.  The  program  is  based  on  a  model  described 
by  Forrest  [Ref.  9:  pp.  35-38].  The  model  requires 
submarine  displacement  as  an  input.  Table  IV-6  displays 
results  when  the  submarine  displacement  is  doubled  in  each 
succeeding  entry. 


Table  IV-6.  Slant  Detection  Ranges  in  Meters  for 
Different  Submarine  Tonnages. 


Displacement 

Signal 

Slant  Detection  Ranges 

in  Meters 

in  tons 

Magnitude 

Cross- 

Square 

OPTEVFOR 

in  oersted 
cm3 

Correlation 

Law 

1000 

1.59x10® 

590 

463 

236 

2000 

3.17x10° 

732 

575 

297 

4000 

6.35x10® 

907 

714 

376 

8000 

1.27x10’ 

1127 

885 

472 

16000 

2.54x10’ 

1402 

1099 

597 

32000 

5 . 08x1 09 

1724 

1363 

753 

As  can  be  seen  from  column  two  in  Table  IV-6,  the  dipole 
moment  is  proportional  to  the  displacement.  Since  the  three 
detection  models  give  a  slant  detection  range  that  is 
proportional  to  the  cube  root  of  the  dipole  moment,  doubling 
the  submarine  displacement  should  multiply  the  slant 
detection  range  by  21^3  (1.26).  This  is  confirmed  by 
comparing  the  slant  detection  ranges  between  the  entries  in 
Table  IV-6.  Doubling  the  displacement  multiplies  the  slant 
detection  range  by  1.24  for  the  cross-correlation  and  square 


law  detection  models  and,  as  expected,  by  1.26  for  the 
OPTEVFOR  detection  model. 


Table  IV-7  lists  the  displacement  in  tons  of  selected 
Soviet  submarines.  The  values  were  taken  from  Combat  Fleets 


of  the  World  1982/1  983  [Ref.  14:  pp.  602  -61  4  ],  This  table 


Table  IV-7 .  Selected  Soviet  Submarine  Displacements. 
Class  Displacement  in  Tons 


Typhoon 
Delta  III 
Yankee 
Echo  II 
Victor  I 
Charlie  I 
Tango 
Foxtrot 
Whiskey 


25-30,000 
10,500-13,250 
8,000-9,600 
5,000-6,000 
4,300-5,100 
4,000-4,900 
3,000-3,700 
1 ,950-2,400 
1 ,080-1,450 


is  presented  solely  for  the  purpose  of  the  information  it 
contains.  The  submarine  magnetic  dipole  moment  program 
should  not  be  expected  to  give  accurate  estimates  of  these 
submarine's  Induced  magnetic  moments,  since  the  program  uses 
a  value  that  relates  displacement  to  magnetic  moment  that  is 
based  on  submarines  of  World  War  II. 


E.  VERTICAL  SEPARATION 

Figure  4.10  shows  three  lateral  range  curves  for  a 
vertical  separation  of  500  meters.  The  OPTEVFOR  detection 
model  lateral  range  curve  shows  only  a  slight  detection 
probability  even  when  the  aircraft  passes  directly  over  the 
submarine.  The  cross-correlation  and  square  law  detection 
model  lateral  range  curves  show  an  Increase  in  lateral 
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law  detection  models  lateral  detection  ranges  increase  with 
an  increase  in  vertical  separation  until  about  500  meters. 


Table  IV-8.  Lateral  Detection  Ranges  for  Different  Vertical 
Separations. 

Vertical  Separation  Lateral  Detection  Range  in  Meters 


in  meters 

Cross- 

Correlation 

Square  Law 

OPTEVFOR 

1  00 

304 

614 

360 

200 

885 

685 

31  8 

300 

942 

720 

23 

400 

974 

724 

500 

980 

699 

600 

974 

629 

700 

936 

262 

No  longer  attains  a  probability  of 
percent. 

detection 

equal  to  50 

A  factor  related  to 

vertical  separation  is  the  effect  of 

ocean  wave  noise  on  a 

MAD  system. 

As  the  altitude  of  a 

magnetometer  is  decreased,  the  magnitude  of  the  ocean  wave 
noise  increases.  Because  of  the  rate  of  this  increase,  for 
a  given  submarine  and  submarine  depth  there  is  a  minimum 
altitude  at  which  an  aircraft  should  prosecute  a  submarine 
using  MAD.  Further  investigation  using  an  ocean  wave  noise 
model  might  be  valuable. 
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V.  CONCLUSIONS 


This  thesis  has  presented  a  comparison  of  three  MAD 
detection  models.  The  cross -corre la tion  detection  model, 
which  models  an  optimum  detector  under  the  conditions  of  the 
detection  model,  yields  the  maximum  detection  range  tor  a 
set  of  given  conditions.  The  square  law  detection  model 
does  not  describe  an  optimum  detector  under  the  conditions 
of  the  model  and  yields  shorter  detection  ranges.  In  the 
stochastic  mode,  with  an  appropriate  choice  for  the 
parameter  AL  that  determines  the  standard  deviation,  the 
lateral  range  curves  for  the  OPTEVFOR  detection  model  become 
similar  to  the  other  two  detection  models.  Detection  ranges 
for  the  OPTEVFOR  detection  model  depend  on  the  choice  for 
the  Operator  Recognition  Factor  (ORF).  With  a  value  of  3 
for  the  ORF,  it  yields  the  shortest  detection  ranges. 
Adjusting  the  magnetic  noise  level  by  an  amount  proportional 
to  the  effective  ORF,  the  cross-correlation  and  square  law 
models  can  be  used  to  describe  the  performance  of  an 
operator  or  an  automatic  detection  system. 

The  magnetic  signal  and  magnetic  signal  plus  noise 
traces  appear  to  have  the  characteristics  of  actual  signal 


traces.  This  suggests  that  the  signal  and  noise  models, 
which  are  the  basis  for  the  cross -correlation  and  square  law 
detection  models,  might  be  useful  for  training  purposes. 


m 


Variations  on  a  set  of  base  case  parameters  were  used  to 
show  relative  changes  in  the  detection  models.  The 
parameters  included:  magnetic  noise,  submarine  and  aircraft 
magnetic  headings,  submarine  displacement,  and  vertical 
separation.  Significant  results  were  the  large  asymmetry  of 
the  lateral  range  curves  under  certain  conditions  and  the 
variation  of  the  magnetic  signal  as  shown  by  the  changes  in 
vertical  separation. 

The  FORTRAN  and  BASIC  programs,  along  with  an  input 
parameter  discussion,  are  included  to  facilitate  the  use  of 
the  three  MAD  detection  models  as  they  are  implemented  by 
the  programs. 
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to  dim  a<2a0),Diaaa),D2aa(a5,K<io0).x0<i00) 

IS  DIM  DSC  100) 

20  DEO 

30  PRINT  “MAX  FREQ") 

40  INPUT  FI 

30  PRINT  i PRINT  "MAX  FREQ  -  "1F1 
SO  Tl-l/Fl 

70  PRINT  "INTERVAL  TIME  -  " I  *  INPUT  T? 

90  G-T7/2/T1 
90  H-INT<Q> 

100  H-H+INT<2t*«:a-H)) 

110  M«2*H+1 

120  IF  M) 200  THEN  70 

130  PRINT  "INT  TIME  -  "1T7 

140  T7-T1*M 

130  PRINT  "ADJ  INT  TIME  ■  "1T7 
1SB  PRINT  "SAMPLE  SIZE  ■  "IM 
170  PRINT  i PRINT  “P/A  RATE  "I 
ISO  INPUT  F2 

ISO  PRINT  « PRINT  "F/A  RATE  «  "lF2 
200  pi»F2*<m-i>*ti/3S00 
210  PRINT  "PP  ■  “1P1 

220  PRINT  l PRINT  "INPUT  DIP  ANCLE  Cl-VES.  0-NO) " II  INPUT  A 

230  IF  A-0  THEN  200 

240  PRINT  "DIP  ANdLE  PHI  “1» INPUT  F 

230  OOTO  420 

200  DEO  IL1-7HL2-100 

270  PRINT  » PRINT  "LATITUDE  "1 

2*0  INPUT  L 

290  PRINT  IPRINT  "LONQITUDE  "» 

300  INPUT  0 

310  PRINT  IPRINT  "LAT  «  "1L«PRINT  "LON  -  " 10 

320  p-siN(o-c2)*GDa<u  «a*caa<o-L2>*caa<u  ih-sincd 

330  U-aiV-H 
340  G08UI  1900 

300  J«J-<90-Ll) «0*K*aiN<J) iH«K*C08U> 

3B0  U-OIV-F 

379  aosua  1900 

380  F-KlO-0lR—<Caa(;Ll>t*0IN<J))ia— (CQBCLD-Hzasu)) 

310  U-Hiv-P 

400  Q0SUU  1900 

410  P-ATN<2*<8IN<J)/Caa<J)>> 

420  PRINT  "PHI  ■  "IF 

430  PRINT  IPRXMT  "DIPOLE  COURSE  "MINPUT  CUPRINT  "DIPOLE  SPEED  "MINPUT  VI 

440  PRINT  IPRINT  "SENSOR  COURSE  "MINPUT  C2IPRINT  "SENSOR  SPEED  "MINPUT  V2 

480  PRINT  "DIPOLE  COURSE  ■  "tCHPRXNT  "DIPOLE  SPEED  -  “»V1 

480  PRINT  "SENSOR  COURSE  •  'MCSiPRXNT  "SENSOR  SPEED  *  "1V2 

470  Ml»V2*SIN<C2)-VtftIN<CJ)lN2-V2*Ca8CC2>-Vl*C0i<Cl) 

480  U*W1 lV"W£ 

490  aOBUS  1900 
800  C0*JiW0"K 
310  D3-W0*T 1*4.  83/9 

320  PRINT  "REL  COURSE  -  "IC01 PRINT  "REL  SPEED  -  "IW0 

830  PRINT  IPRINT  "INPUT  DIPOLE  MOMENT  a«YES.  0-NO) " M INPUT  AA 

840  IF  AA-0  THEN  300 

830  PRINT  IPRINT  "MAQNITUDE  P  "MINPUT  PlPRINT  “HOR  ANGLE  w  "MINPUT  A 
3B0  PRINT  IPRINT  "VERT  ANGLE  OMEGA  "MINPUT  B 
870  OOTO  340 

800  PRINT  IPRINT  "INPUT  EARTH  FIELD  <  1-YBS,  B-NO)  "MINPUT  AA 
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590  IF  AA-9  THEN  620 

600  PRINT  « PRINT  "EARTH  FIELD  "I i INPUT  El 
610  GOTO  630 

620  & 1 -70000/ 9QR ( 3-COS <  F  )  -COS  <  F!) + 1  > 

630  PRINT  » PRINT  "EARTH  FIELD  -  'ME! 

640  (14-01 113-0  mS-0 

630  PRINT  iPRINT  “INPUT  PERM  MOMENTS  <1 -YES. 0-NO)  “it  INPUT  AA 
6E0  IF  AA-0  THEN  630 

670  PRINT  IPRINT  “LQNQ  MOMENT  HU INPUT  M4» PRINT  “TRAN  MOMENT  "I  I  INPUT  MS 
SS0  PRINT  "VERT  MOMENT  “li INPUT  M« 

630  PRINT  iPRINT  “LONQ  MOMENT  -  " 1M4I PRINT  "TRAN  MOMENT  -  "iMSiPRINT  “VERT  MOMEN 
T  -  “IMS 

700  Kl-7, 3IK2-1. 6IK3-1. B 

719  PRINT  IPRINT  '-DISPLACEMENT  " I 

/20  INPUT  N1 

730  PRINT  iPRINT  “DISPLACEMENT  -  "INI 
733  NN1-N1 

740  M9-E1-K3*N1-SIN(F) IM3-M9+M6 

750  Mfl-Et*CaS(F)»Nl“'<Kl#C0Si:Cn*caB<Cl)+K2*SINt:Cl)*SIN(Cl)) 

760  Ml-M4*SIN<Ci >-MS*CQS  CCt ) IM2-M4-C0S (Cl ) -MS-SINCCl ) 

770  M7"Si*<CaS<F>>HKl-tf2>-Nl*9IN<Cl>*C03<:Cl)> 

700  M 1 -M7+M X I M2-MS+M2 
790  U-MUV-M2 
300  QOSUB  1300 
310  A-JIU-M3IV-K 
920  QOSUB  1300 
330  P-KIB-J 

340  PRINT  IPRINT  “P  -  “IPI PRINT  “W  •  “lAl PRINT  “OMEGA  -  “IB 

350  X-Pl 

963  aaSUB  1320 

670  VS-Y 

900  V7-M*<1-2/9/M4Y*0QR<2/9/M>>a3 

390  PRINT  iPRINT  “VERT  SEPARATION  "»* INPUT  l 

900  PRINT  IPRINT  "VERT  SEPARATION  ■  "1Z 

910  PRINT  iPRINT  “NOISE  “1 

920  INPUT  81 

923  92-01-4 

930  PRINT  iPRINT  “NOISE  -  *1111 

940  PRINT  IPRINT  "MAX  LATERAL  RANOE  “M INPUT  R0 

950  PRINT  iPRINT  “NUM1ER  OF  INCREMENTS  “H INPUT  N7 

9C0  PRINT  iPRINT  “MAX  LATERAL  RANOE  «  “1R8«PRINT  “NUMBER  OF  INCREMENTS  -  “1N7 

9E4  PRINT  iPRINT  “  ORF  "IHNI'UT  ORF 

9EE  PRINT  iPRINT  “  ORF  •  “10«F 

970  D4-RE/N7IN0-2-N7 

174  PRINT  IPRINT  “  ALPHA  "M  INPUT  AL 

978  PRINT  IPRINT  “  ALPHA  •  “1AL 

910  LS— R0 

910  FOR  E-0  TO  N0 

1000  X0-L9 » X0  <E) -LS 

1010  aaius  1B29 

1020  LS-LS4D4 
1030,  QQIU0  1,800  , 

UMwwr  a 

1044  aOTO  1230 

1050  ORAPHICS  3 i COLOR  1 

1000  XX-INT1310/N0) 

1005  X0<0>-0 

1070  FOR  1-0  TO  N0 

1000  DlCD-INTat-DiaJJ-iSO) 

1010  D2a)-INT<<:-D2a>>H<180) 

1093  03C I)-INTt (i“D3< I) )-lS0) 

1190  X0<I-fO-X0<I)+XX 
1110  NEXT  I 
1120  PLOT  X0<3> i 01 <05 
1130  FOR  1-1  TU  Ni 
1140  DRAWTO  XOO.DKI) 


'vIvvTjv^vVVj f. 
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USB  NEXT  I 

116(9  PLOT  X0(05.D2(0) 

1170  FOR  1-1  TO  N0 
1160  DRAWTO  XOCD.QSa) 

1190  NEXT  I 

119?.  PLOT  X0(05  t  D3<05 

1194  FOR  !•!  TO  N8 

1195  DRRWTO  X0<I>,O3<I) 

1190  NEXT  I 

1200  PRINT  "PD  FOR  X  FROM  "*-RBl"  TO  "IRS 
12S0  GOTO  1810 

1290  PRINT  "FOR  HARD  COPY  ENTER  •  1"  "l  i INPUT  CO 

1233  IF  CCOl  THEN  GOTO  1030 

1300  LPRINT  "MAX  FREQ  -  'MFl 

1310  LPRINT  "AOJ  INT  TIME  -  *1T7 

1320  LPRINT  "SAMPLE  SIZE  -  'MM 

1330  LPRINT  "F/A  RATE  -  "*F2 

1333  LPRINT  "PF  •  "1P1 

1340  LPRINT  "LAT  •  ",1L 

1330  LPRINT  "LON  -  'MO 

1380  LPRINT  "PHI  «  'MF 

1370  LPRINT  "DIPOLE  COURSE  -  "IC1 

1380  LPRINT  "DIPOLE  3 PEED  -  “1V1 

1390  LPRINT  "SENSOR  COURSE  -  "*C2 

1400  LPRINT  "SENSOR  SPEED  ■  "*V2 

1410  LPRINT  "REL  COURSE  -  " ICO 

1420  LPRINT  “REL  SPEED  -"1W0 

1430  LPRINT  “BARTH  FIELD  -  "Ittl 

1440  LPRINT  "LONO  MOMENT  «  »tM4 

1430  LPRINT  "TRAN  MOMENT  -  "IMS 

1480  LPRINT  "VERT  MOMENT  «  «1M*i 

1470  LPRINT  "DISPLACEMENT  -  "!NN1 

14S0  LPRINT  "P  -  "IP 

1490  LPRINT  "W  «  "*A 

1300  LPRINT  "OMEOA  -  "IB 

1303  LPRINT  “VERT  SEPARATION  -  "1Z 

1310  LPRINT  "NOISE  -  "101 

1813  LPRINT  "MAX  LATERAL  RANGE  -  MRS 

1320  LPRINT  "NUMBER  OF  INCREMENTS  *  "IN7 

1323  LPRINT  • LPRINT  "LTN  RNQ  PD(CC)  P0<SL)  PD(OPT)" 

1330  LS— SB 
1833  FOR  1-0  TO  N« 

1340  LPRINT  LSI"  "101(1)1*  "102(15*“  "»D3(I5 

1343  L3-L3+04 

1330  NEXT  I 

1360  QOTQ  1030 

1803  PRINT  "END" 

1810  END 
1620  U-X0IV-Z 

1630  aaauii  isoa 

1640  D-JIH0-K 

1642  RH-<0. 1*P/ <0RF*62> ) *0. 333 

1643  8XQ-AL*RH 
1843  X-(RH-H0)/'6ia 

■nraraoainniiB— 

1647  DSCEi-V 

1630  S0«caacl5*ca8<c0-H> »JB"CaB<05*Caa<B>*BIN<C0-A>-SIN<D)*9IN<B> 

1860  N0»-<81N<D5#CQS<1)|MHN<C0-A)  )-C0S<D5|"SIN<B> 

1870  »l-Caa<F)*COS<C05iJl-Cas<D5*Ca3<F5*8IN<C05~SIN<D>*6,IN<F5 
1660  Nl**-<aiN<D>*C08(F)*«IN<C(’)>  >-CQS<D)*8IN<F5 
1690  K1»P/10/H0*3 

1700  A2-2*B0*Bl-J0*JlH'»0#NlkAl«3*<N0*8l-t-B0*Nl>  i  A0»'2*N0*N1'-B0*B1-J0*J1 
1710  90-0 

1720  FOR  1-0  TO  M-t 
1730  9"m;I-i:m-1)/25*03 
1740  Q-S/HB 


nan  G-im  ;'Q*q>a2.  s 
17SG  g«<A2*q*q+ai*q+a0>*g 
1772  G<i:>-G!S0-S3+G*a 
1780  NEXT  I 
1730  RETURN 
1800  K2*SQR<a0) 

1310  K<E)«K2 

1820  V8»~V6fKl’*,8QR<S0)  /S1 

1330  LO-K1*K1*S0/ <S1*>811 »A3-M+L0»B3-1+L0/ (M+L05 
1840  V8»-8QR<2*V7/B3>+SQR(2*A3/B3-1) 

18S0  X-V8IGGSUB  1990 
1860  01<E)-Y 
1370  X-V9IG0SUB  1390 
1380  D2<E>-Y 
1830  RETURN 

1800  K»3QR<U*U+V*V)«IF  K«0  THEN  J"0« RETURN 
1908  UK-U/KIVK-W/K 

1907  IF  UX> 0.999999  AND  VK>0. 9S9399  THEN  J«0i RETURN 

1908  IF  UK>0, 339399  THEN  J— ATN<VK/9QR<-VK*VK+i > )+30« RETURN 

1909  IF  VK>  0.999999  THEN  J-01 RETURN 

191.0  MM-ATN<UK,'SQR<-UK^UK+1))'.J— ATN<VK/3QR<-YK**.VK+1))>90«IF  MM<0  THEN  J-360-J 
<913  RETURN 

1920  Y»X  5  IF  5O0.3  THEN  Y-i-Y 
1930  Y»8GR  <  LOG  < 1 / V/  V  >  ) 

1940  Uft-2. 3158 1 7 1 Gl-0. SC2853 1 G2-0.  0 1 0323 

1930  Hl-1.  432788KH2-0. 183269«H3-1. 308E-03 

1860  Y-Y~  < G0+Y*  <G  J.  +G£*Y  )  )  /  <1 + Y*  <H1  Y*  <  H2+H3*Y )  >  > 

1370  IF  X)  0. 3  THEN  Y— Y 
1980  RETURN 

199B  Y-XUF  X<0  THEN  Y— Y 
2000  V)*l/a+0.2316419t-Y) 

2010  Qi«C.  31938133«Q2— 0.  336363782 »Q3*1. 78 1 47793 « Q4» “ i. B2 12233 3 «Q5-1. 33027442 
2020  IF  Y>  24. 23  THEN  YMStGOTO  2070 
2023  Pl-3.  14139263 

2030  Y«SXr<-<Y*Y/25 ) /SQR  <  2*P I) <  Q1 +W*  <  Q2+U*  <  Q3+W*  < 04*14*03 ) ) ) ) 

2070  IF  X>0  THEN  Y-l-Y 
2073  Y-<INT<1000«*Y) 1/10000 
2080  RETURN 


P 
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MAX  FREQ  -  1.13 

ADJ  I  NT  TIME  *>  20. 33333335 

SAMPLE  SIZE  -  37 

F/A  RATE  -  3 

pr  -  9.016666666G 

CAT  -  30 

LON  -  60 

PHI  «  39.40076979 
DIPOLE  COURSE  «  43 
U I POLE  SPEED  -  10 
SENSOR  COURSE  -  31S 
SENSOR  SPEED  -  220 
REL  COURSE  -  312.397439 
REL  SPEED  -220.227133 
EARTH  FIELD  -  32306.3313 
LPNO  MOMENT  -  0 
TfirtN  MOMENT  -  0 
VERT  MOMENT  -  0 
DISPLACEMENT  -  4000 
P  -  624S94392 
m  -  32.  S37307S1 
OMEQA  -  27. 11173179 
VERT  SEPARATION  -  200 
NOISE  -  0.1 

MAX  LATERAL  RANGE  «  1330 
NUMBER  OF  INCREMENTS  -  13 


* 


LTP  RNO 

PD(CC) 

POOL) 

1 

-1300 

0. 0344 

0.0138 

0 

-1400 

0. 0696 

0.0166 

0 

-1300 

0. 0948 

0.0181 

0 

-1200 

0. 1391 

0.0209 

0 

-1100 

0.2212 

0.0266 

0 

-ia'00 

8. 3763 

0.0403 

0 

-900 

0. 6399 

0.0801 

0 

-800 

0.9191 

0.  2233 

0 

-700 

0.  999 

0. 6743 

0 

-600 

1  0. 

9973  2E-04 

-300 

1  1 

0.0135 

-400 

1  1 

0.1603 

-300 

1  1 

0. 367 

-200 

1  1 

0. 8868 

-100 

1  1 

0. 9774 

0  1 

1 

0. 3898 

100 

1  1 

0. 9774 

280 

1  1 

0. 3868 

300 

1  1 

0.  367 

400 

1  1 

0. 1603 

300 

1  0.99S9  0.0133 

600 

0. 9997 

0. 7C34 

2E- 

700 

0.9137 

0.2189 

0 

800 

0.3773 

0. 2677 

0 

900 

0. 3014 

0. 0331 

0 

1000 

0. 1644 

0. 0226 

0 

1100 

0.  1003 

0.0184 

0 

1200 

0. 0683 

0.0163 

0 

1300 

0. 0308 

0. 0136 

0 

1400 

0. 0403 

0.0131 

0 

1300 

0.  034 

0.0143 

0 

PD (OPT) 


l .  t 
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